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Many correlated materials display a quantum critical point between a paramagnetic and a spin- 
density wave (SDW) state. The SDW wave vector connects points, so-called hot spots, on opposite 
sides of the Fermi surface. The Fermi velocities at these pairs of points are in general not parallel. 
Here we consider the case where pairs of hot spots coalesce, and the wave vector (tt, tt) of the SDW 
connects hot spots with parallel Fermi velocities. Using the specific example of electron-doped 
cuprates, we first show that Kanamori screening and generic features of the Lindhard function 
make this case experimentally relevant. The temperature dependence of the correlation length, 
the spin susceptibility and the self-energy at the hot spots are found using the Two-Particle-Self- 
Consistent theory and specific numerical examples worked out for band and interaction parameters 
characteristic of the electron-doped cuprates. While the curvature of the Fermi surface at the 
hot spots leads to deviations from perfect nesting, the pseudo-nesting conditions lead to drastic 
modifications to the temperature dependence of these physical observables: Neglecting logarithmic 
corrections, the correlation length ^ scales like 1/T, namely z — 1 instead of the naive z — 2^ the 
(7r,7r) static spin susceptibility x hke 1/a/T, and the imaginary part of the self-energy at the hot 
spots like T^/^. The correction T^^ ^ rp2,/2 ^j^^ Korringa NMR relaxation rate is subdominant. 
We also consider this problem at zero temperature, or for frequencies larger than temperature, using 
a field-theoretical model of gapless collective bosonic modes (SDW fluctuations) interacting with 
fermions. The imaginary part of the retarded fermionic self-energy close to the hot spots scales as 
— cj^/^ logcj. This is less singular than earlier predictions of the form —uj logcj. The difference arises 
from the effects of umklapp terms that were not included in previous studies. 

I. INTRODUCTION 

Quantum phase transitions between a Fermi liquid and magnetic phases have been a subject of experimental and 
theoretical investigations for several decades . The transition to a spin-density- wave (SDW) in particular is relevant 
to problems of current interest. In the cuprates, Daou et al. argued that the Fermi surface change associated with 
this transition is a key to understanding anomalous normal state properties. Recent studies in the pnictides , heavy 
fermion materials and organic superconductors focus on the relation between the SDW and superconductivity. 
In fact, strong experimental similarities between quantum critical behavior in the organics, pnictides and cuprates 
have been pointed out recently . 

The electron-doped cuprates have provided an early example where quantum critical behavior has been inferred 
from the temperature dependence of resistivity at low temperature. It was measured to be linear from 35mK to 
lOK in Vi2-xCexCn04^-s (PCCO) at doping x = 0.17. More recent transport ^ and thermopower measurements 
also suggest the presence of a quantum critical point at a similar doping. The precise nature of the quantum critical 
point remains however unclear, as thoroughly discussed in Ref. . For example, it has also been suggested that the 
quantum critical point coincides with the onset of superconductivity in the overdoped regime . 

Here we study quantum critical behavior associated with the transition between a SDW phase and a Fermi liquid 
when the wave vector (tt, tt) of the SDW connects hot spots with parallel Fermi velocities. The two Fermi surfaces 
connected by (7r,7r) in this case are tangent to each other, as shown in Fig. 1. On a spherical Fermi surface, the 
Fermi wave vector would satisfy the condition 2kF = (7r,7r). We call this pseudo-nesting . We will explain why the 
QCP can be located at, or close to, the pseudo- nesting filling 71^ This occurs naturally in the one-band Hubbard 
model for the electron-doped cuprates and we will perform some of our calculations specifically for this case, although 
the frequency and temperature dependencies that we find are valid more generally. The methods that we describe 
below can be applied to electron-doped cuprates because these materials are described by a Hubbard model in an 
intermediate coupling regime where one can neglect effects induced by Mott physics . 

The theory of Hertz ' ' and Millis has formed the basis for much of the work on quantum critical phenomena. 
In this approach, fermions are integrated out and an effective bosonic theory for the collective modes is studied using 
standard renormalization group methods that can be taken to high order^''. It has been pointed out by Abanov 
and Chubukov^'^ that for a commensurate SDW at the upper critical dimension, namely (i = 2 for 2; = 2, all the 
coefficients of the bosonic theory are singular so that one must treat simultaneously the bosonic collective modes and 
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FIG. 1: The dashed hnes indicate the magnetic Brillouin zone with ordering vector Q = (tt^tt). The arrows gives examples of 
pseudo-nesting conditions, namely of points such that 2kF is equal to the antiferromagnetic wave vector. The (tv^tv) ordering 
wave- vector is defined with respect to the (qx, Qy) coordinate system. In the field theory approach introduced later, we work in 
the rotated (kx^ky) coordinate system. 



the fermions. Metlitski and Sachdev^^ have reexamined this problem and obtained the non-Fermi liquid behavior at 
the hot spots, and shown that the bosonic SDW spectrum does not obey dynamic scaling with z = 2 but instead that 
a super-power-law form is obtained. They have also thoroughly discussed the failure of the expansion at higher 
order, leading to a strong-coupling problem. More recently, it has been argued that non-Fermi liquid corrections 
are also important away from the hot spots. 

An alternate approach is the self-consistent-renormalized theory of Moriya^^'^^. It is in the universality class of the 
spherical model and as such its critical behavior will not be exactly that expected for the 0(3) model. Nevertheless, 
it can be accurate away from the critical point and provide leading order estimates for the exponents. The Two- 
Particle- Self- Consistent (TPSC) theory^^''^"^ is a related approach that has critical behavior similar to that of Moriya, 
including logarithmic corrections . It has the advantage that although it is an approximate solution to the Hubbard 
model, it is quantitatively very close to benchmark Quantum Monte Carlo results . It is non-perturbative, does not 
include any phenomenological parameters, has internal consistency checks, and satisfies a number of exact results. 

Previous theoretical studies have mostly been done for the case where the quantum critical point occurs when the 
Fermi velocities at the hot spots that are connected by the SDW are not parallel, unlike the case of parallel velocities 
we consider here (see Fig. 1). Such a case of parallel Fermi velocities is generic in one dimension, but at first sight 
appears as an accident in two dimensions, because upon translation by the (tt, tt) SDW wave vector, the Fermi surfaces 
touch at only one point. If the surfaces were flat, we would recover the case of perfect nesting encountered in one 
dimension. The curvature here provides a cutoff, and so we refer to the situation with parallel Fermi velocities as 
"pseudo-nesting"^^. Such Fermi surfaces have also been studied in three dimensions'^, but in the presence of this 
pseudo-nesting the spin susceptibility is singular in two dimensions and the analysis must be redone. Altshuler, 
loffe and Millis first looked at the case where the instability is at 2kF, hence connects parallel segments of the Fermi 
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surface, but the SDW wavelength is not commensurate with the lattice. They found that the transition is weakly 
first order with an intermediate scaling regime when the SDW wavelength is close to (7r,7r). The scaling regime was 
obtained in a systematic expansion in a number proportional to the inverse number of fermion flavors. Krotkov and 
Chubukov ' found different results for the self-energy. Here we consider only the commensurate case. Some of our 
results differ from those of previous authors because they overlooked the significance of the umklapp process shown 
by the top double arrow in Fig. 1. 

We use two different approaches. We obtain finite temperature results appropriate for the SDW quantum critical 
point of electron-doped cuprates using TPSC. Then a field-theory for the spin-fermion model allows us to find the 
finite-frequency zero temperature results and some finite temperature results. The results of both approaches are 
consistent. 

The rest of this paper is organized as follows. In Sec. II we present the model along with general arguments 
suggesting why one should expect the antiferromagnetic quantum critical point to be located close to pseudo-nesting, 
namely at a filling where the (tt, tt) wave vector connects parts of the Fermi surface that are tangent, or equivalent ly 
with parallel Fermi velocities. Sec. Ill contains the finite temperature results. They are obtained with TPSC, which 
is described in Sec. Ill A. Analytical results for the behavior at the QCP are illustrated with numerical examples 
appropriate for electron-doped cuprates in the subsections of Sec. IIIB. First the low temperature asymptotic form 
of the susceptibility is obtained in Sec. IIIB 1, then the critical behavior is obtained in Sec. IIIC. Zero temperature 
finite-frequency results are treated in Sec. IV with field theoretical methods. The Lagrangian appears in Sec. IV A 
followed by sections on the polarization bubble (spin susceptibility) IV B, on the electron self-energy IV C and on 
the irrelevance of the quartic term IV D. An appendix on vertex corrections A appears after the summary in Sec.V. 
Consistency between TPSC and field theory results are pointed out in the field theory Sec. IV. 



II. MODEL AND QCP FOR ELECTRON-DOPED CUPRATES 

In this section, we introduce the model and give generic arguments why we expect the quantum critical point to 
often be located close to the filling where translation by the antiferromagnetic wave vector leads to Fermi surfaces 
that are tangent to each other, as illustrated in Fig. 1. 

While the frequency and temperature dependencies that we find do not depend on details of the model, specific 
numerical examples at finite temperature calculations will be performed on the two-dimensional t — t' — U Hubbard 
model on the square lattice at weak to intermediate coupling. The model is given by 

H = -^ tij{cl^Cj^a ^ h.c.) ^U^fii^^rii^l (1) 

where tij are the hopping integrals, z, j are the site index, a is the spin label, ^ and q^^t are the particle creation 

and annihilation operators. Doubly occupied sites cost an energy U and n^^^r = cj^Q^cr- Units are such that h = 1, 
ks = I and lattice spacing is unity. The kinetic energy of a single-particle excitation in momentum space is obtained 
from 

£k=|^-E«''-^-'"^^'*^.j-M« (2) 

with the sum over j running over all neighbors of any of the sites i. The chemical potential /i^^^ is chosen so that we 
have the required density. 

One can explain on general grounds the filling where the QCP is likely to occur. Figure 2 displays the Lindhard 
function ^g^^^ along the Qy direction for different fillings. Its maximum is at (7r,7r) so it is symmetric in and Qy. 
There are two remarkable features. First, below a certain doping nc, the maximum value is almost independent of 
filling'^^'^^, and second it falls rapidly as soon as the filling exceeds ric. The filling ric corresponds to the point where 
the Fermi surfaces joined by (tt, tt) touch instead of intersecting. A 0.3% change in filling leads to almost 10% drop in 
value of the susceptibility. If we consider a simple Stoner criterion for the transition, we would conclude that if U takes 
the value Uc = ^ 2.6, then the QCP would be close to this filling ric. This does not require fine tuning because 

the value of U that should enter the Stoner criterion is the value renormalized by Kanamori-Briickner screening 
This renormalized value becomes essentially U independent when U becomes of the order of the bandwidth because 
the two-body wave function creates a cusp to minimize double-occupancy ' and the renormalized interaction cannot 
become larger. This maximum renormalized value in TPSC, Ugp^ takes a value near Uc ^ 2.6. In addition, in TPSC 
the value of Usp self-consistently adjusts itself to the value necessary to prevent a finite temperature phase transition 
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FIG. 2: Lindhard function near the 2k f point as a function of doping for U = 6, t' = —0.175 and = 0.05, values that are 
appropriate for electron-doped cuprates. The rapid fall with filling larger than 1.201, close to the critical filling, is apparent. 

on the SDW side of the QCP. Although, at sufficiently low temperature, details will start to matter and one needs to 
start to tune the value of U to find the QCP precisely at nc, there is an intermediate temperature scale that can be 
quite broad where fine tuning is unnecessary. 

III. FINITE TEMPERATURE RESULTS AND TPSC 

In this section, we use the non- perturb at ive Two-Particle Self-Consistent (TPSC) approach'^'^''^^. This approach 
respects the Pauli principle, the Mermin- Wagner theorem and conversation laws. It also contains quantum fluctuations 
in crossed channels that lead to Kanamori-Briickner screening. It is valid in the weak to intermediate coupling 
regime {U ^ 6t) and not too deep in the renormalized classical regime where a pseudogap is observed. It has been 
benchmarked on Quantum Monte Carlo calculations on the Hubbard model. 

TPSC has been shown to be in the N = oo universality class \ It has the same critical behavior as Moriya theory 
and hence has the same logarithmic corrections . These logarithms have the same functional form as those of the 
renormalization group asymptotically close to the quantum critical point, but in TPSC and in Moriya theory the 
mode-mode coupling term does not flow, hence the corrections may differ in the details from the renormalization 
group ^ Quantum critical behavior of the susceptibility and of the self-energy in the closely related spin-fermion 
model has been discussed by Abanov et al.^^^. 

It has been argued from detailed comparisons of numerical calculations with experiment^^'^^''^^'^^''^^ that strong- 
coupling physics is not important for electron-doped cuprates, at least not too close to half-filling. Hence, TPSC is 
appropriate to study these compounds. It gives a satisfactory description of ARPES data^^, and the temperature 
T*, where the pseudogap seen in ARPES opens up experimentally, corresponds to that where the antiferromagnetic 
correlation length coincides with the thermal de Broglie wavelength , as predicted for two-dimensional precursors of 
three-dimensional long-range order . 

Hence, all the numerical results are presented in units where t = l,/c^ = l,fi=l (with z component of spin defined 
by — ni) for values of the Hubbard model hopping parameters appropriate for electron-doped cuprates, namely 
second and third nearest-neighbor hopping = —0.175 and t^^ = 0.05 . Interaction strengths U = 6 and U = 5.56, 
again in the range appropriate for electron-doped cuprates, will be considered. 

We first present the formalism and then give analytical and numerical results for the QCP. 
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A. TPSC 

Given the Hubbard model parameters, TPSC has no adjustable parameter. Irreducible vertices are obtained self- 
consistently and in such a way that the Pauli principle and conservation laws are obeyed. The formal derivation is 
given in Refs. , . Here, we simply present the equations that are solved. 

In TPSC, the retarded spin Xspi^.^^) and charge Xch{^.^^) susceptibilities are written as 

Xch q, ^) = , ^ u, mf — ; ' (^) 

where X^^Hq^^) non-interacting retarded Lindhard function at wave vector q and angular frequency uj 

Here, / (sk) is the Fermi function (e^"^/-^ + l) ^, T is the temperature and N is the total number of sites. The 
effective spin interaction Usp is evaluated without adjustable parameter using the ansatz^^'"^^ 

U (n^n^) = Usp{n^){n^) (6) 

with the local-moment sum rule that follows from the fluctuation-dissipation theorem 
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I . , duj dPq 2 // / N /^N 

n - 2{n,n,) = - ^Y^^X.M-) (7) 

where Xsp ~ ^^Xsp (^t^i) double occupancy. We dropped the site index using translational invariance 

and we used the Pauli principle to write 

52 = ((n^ -71^)2) =n-2(ntn;). (8) 
Similarly the irreducible vertex Uch entering Xch{Q) is found using a sum-rule that is the analog of Eqn.(7) for spin: 

. 2 f°° dio f°° Cpq 2 „ , , 

The crossing symmetric self-energy is obtained from 

ES')(fc) = C/n_, + f ^ E l^UspXspiq) + UehXchiq)] G^^\k + q). (10) 

The superscript (2) reminds us that we are at the second level of approximation. G^a^ is the same Green's function as 
that used to compute the susceptibilities x^^H^)- Charge fluctuations Xch{Q) are included in numerical calculations 
but they are neglected in the analytical results because they are small. 

Since the self-energy is constant at the first level of approximation, this means that G^^ is the non-interacting 
Green's function with the chemical potential that gives the correct filling. This chemical potential /i^^^ is slightly 
different from the one that we must use in (C^^^) = iuj^ + /i*^^^ — — S*^^^ to obtain the same density'^^. 

Unless otherwise specified, all the numerical results below are obtained using the Matsubara frequency version of 
equations (3) to (10) without any approximation, hence they are valid at arbitrary distance from the quantum critical 
point. 



B. Analytical results with numerical examples for electron-doped cuprates 

We begin below with the Ornstein-Zernicke form of the spin susceptibility that is usually valid when the correlation 
length is large. The case where there is perfect nesting leads us naturally to the pseudo- nesting condition relevant for 
this paper. The situation, however, is not as simple as usual since the Ornstein-Zernicke form for the spin susceptibility 
is incorrect in our case, as we will explain. The self-energy is treated at the end of this section. 
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1. Ornstein-Zernicke form for the susceptibility 



a. General case When the correlation length is large, one usually assumes that the denominator of the spin 
susceptibility can be expanded around the wave vectors Qd, where the maxima in x'-^^ occur in d— dimensions. One 
then obtains 



x"p(q,w) 



Vro 



(c^/ro) 



2 ' 



(11) 



where q is measured with respect to the wave vector Q^^ where the spin suseptibihty is maximum ((tt, tt) in our case). 
Defining Umf = {Qdfi) as the value of the interaction at the mean-field SDW transition, the other quantities 

in the previous expression are 



e 

5U 

1 



Ugp 



= u, 



mf 



1 



S^x^l) (q,0) 



2x(i) (q, 0) dq^ 
1 (?x(i)"(q,w) 



Co'x(i) (q,0) 



(12) 
(13) 
(14) 

(15) 



In the expression for the spin susceptibility, the denominators are expanded around the (tt, tt) wave vector. 

To obtain analytical results for the imaginary part of the self-energy Y^^'^^"^ {kp-,^] T = 0) in Eqn.(lO) we use the 
spectral representation for the susceptibilities and for the Green's function, perform the sum of the internal Matsubara 
frequency and then the analytical continuation neglecting the charge fluctuations, to obtain 



where the component of q parallel to the Fermi velocity v^?, is obtained from the solution of the equation 

For all Fermi wave vectors, where ^k^+Qd — ^^e above equation reduces to 

■ q_ uj -\- uj' 

where is the Fermi velocity in the hot region, i.e. where SkiT^+Qa — 0- 

In the asymptotic form of the spin susceptibility Eqn.(ll), the wave vector appears only in the form qy , ^^^^ 
keeping this general form in the equation for the self-energy Eqn.(16), we obtain 



S"«(k^,c.) = -^^ I I ^[n{oj') + f{uj + oj')]x':,{^^,{o^ + o^')/vF;oj',T). (18) 

Normally, one expects to be a temperature independent constant of the order of the lattice spacing and Fq/^o 
to be a constant of the order of the Fermi energy. In the case of perfect nesting, or of pseudo-nesting, this is not the 
case. 

b. Perfect nesting 

Although the case we are interested in does not correspond to perfect nesting, understanding that case first will 
facilitate our task later. There is perfect nesting when the equality ^k = — ^k+Q^ is satisfied for all wave vectors, with 
Qd the nesting wave vector. This case was treated by Virosztek and Ruvalds . The quantities and Fq that are 
usually assumed temperature independent, here become temperature dependent. We show this below. 

For perfect nesting, the Lindhard function becomes 



(19) 
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so that changing to an energy integral we have 

x(i)"(Qd, w) = -TT y dENd (E) (1 - 2/ (E)) d {CO + 2E) (20) 
= .iV,(|)tanh(|;) (21) 
where Nd (E) is the density of states. The real part at zero frequency on the other hand is given by 

In two dimensions, there is a weh known logarithmic divergence of the density of states Nd{^) at the van Hove 
singularity. Neglecting this logarithmic divergence that appears only for a special filling in the hole-doped case, we 
take N(i{^) as a constant. In that case, integrating by part and replacing the upper bound by infinity in the convergent 
integral we are left with 

x(^)(Qd,0) = iVd(0)(^ln(x)tanh(x)|o^-y--^^dx^ « iV, (q) In + B (23) 

where 5 is a temperature independent constant. We also have 

x(i)"(Qrf, w) « 7riVrf(0)tanh (^) . (24) 

These results suggest that the quantity Fq defined by Eqn.(15) scales as 

ro-eo'^ln(A/T). (25) 

Following Ref. , we move on to demonstrate analytically that in Eqn.(14) scales as ~ The 
dependence, fundamentally comes from the second derivative of X^^\Qd^ 0) ^ N (0) In (^) in Eqn.(23). We shall now 
make this argument more rigorous. Keeping for a while a general notation where i is some direction in the Brillouin 
zone, and q is measured with respect to the center of the zone, one can write 

d\^'H^fi)__^f d'^k d'C fde^+^Y_^f d-'k aV+q ^26) 



Bzi'^^Ydel^^ V % ; yB^(27r)^aek+q Bq^ 



where 

/(eic+q) - /(eic) 



Ck+q — Ck 



Measuring q with respect to Q^^ we evaluate the above second derivative at q = 0. As before, for perfect nesting we 
have 

^/ ^ N 2/ (ek+Qd+q) - 1 tanh (ek+Q^+q/2T) 

0(ek+Q^+q,ek) = — = — (27) 

^Ck+Qd+q ^Ck+Qd+q 



1 /^^k+Qd+q\ 



(28) 



The last equation shows that C scales as times a function of ek+Qd+q/^- In the integrals, the derivatives of the 
type will not introduce singular terms in temperature. Hence, replacing them by some average value in the 

Brillouin zone, we can change the integration variable to energy and the most singular terms in temperature will come 
from 



a^xd) (q,0) 



q=0 



■ JdEN.iE)^ {v,f (29) 



^ / dENa {E) ^AJ_r2 (^^)2 (3^^ 
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FIG. 3: The Fermi surface at the critical fiUing Uc touches the antiferromagnetic Brillouin zone. The important integration 
region is over a rectangle of thickness T and width \/T. The critical chemical potential jic where the antiferromagnetic zone 
boundary touches (7r/2,7r/2) on the Fermi surface is the solution of — 2t^^(cos(7r) + cos(7r)) — /ic = 0. This corresponds to a 
filling Tic = 1.2007 for t' = -0.175 and t" = 0.05. For [7 = 6, the critical fihing that we find, ric = 1.20096, is slightly larger. 



Neglecting the energy dependence of the density of states, we are left with 

d\^'^ (q,0) 



q=0 



1 



dxN, (0) (v,)^ 



Using the definition of Eqn.(14), and the result for Fq Eqn.(25) above, we have that 



^0^ 



1 



To 



ln(A/r) 



(31) 



(32) 



c. Pseudo-nesting In this subsection we show that, for the pseudo-nesting case, the main contribution to the 
Lindhard function at Q^^ has the same form as in the perfect nesting case except for a temperature dependant 
prefactor. The previous calculation illustrates that the main contribution to the quantities of interest, and Fq, 
come from nested regions of the Fermi surface. In the pseudo-nesting case ' illustrated in Fig. 3, the Fermi surface 
displaced by Q^^ just touches the original Fermi surface, with the Fermi velocities of the two surfaces that are parallel 
at the touching point. Similarly to the perfect nesting case, the most important contribution to the integral for the 
Lindhard function around Q^^ comes from the regions in /c-space connected by Q^^, with a width around the Fermi 
surface that corresponds to an energy range T. Now, imagine that we divide the integral over k near one of 

those points of the Fermi surface, for example = (f , f ), into two components, k^ parallel to the Fermi velocity, 
i.e. perpendicular to the Fermi surface, and k|| parallel to the Fermi surface at k^. For k near that region of the 
Fermi surface, we can write in two dimensions 



^k||+k^ 



^gk||+k^ 

Skx 



1 9^s:k||+k 



-Ski 



where we have measured wave vectors with respect to k^ and used the fact that Vkn^kn+kj 
measures the curvature of the Fermi surface. From that approximation we have Sku+k^+Qd — 



thus 



(33) 



= 0. The quantity n 



^k||+k^ — ^k||+k^+Qd — "^vpSkj^ . (34) 

Since the terms in cancel out in that expression, this approximation is valid if the next term in the series of 
^k||+k^ is negligible compared with the first one, namely if 



1 a^sk, 



dkl 



5kl 



(35) 
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and since VF5k± < T, we have the fohowing upper bound for the temperature, 



1 ^^^^\\+^± 

2 



(36) 



Now for the Fermi function, over a region around k^, we have 



/(£k+Qj -f{-VFSkj_) 



dfie) 



ds 
dfie) 



e—VF5k± 



ds 



e=—Vf6k_\_ 



(37) 



Since the derivative of the Fermi function is even in 6, 

/(£k)-/(ek+Qj -2/(^;f,5fci)-l. 
The region where this is vahd is given by the condition 



2 ds^ 



£—VF5k± 



which is satisfied if Kdk'^^ ^ T. Therefore, from (34) and (38), the Lindhard function takes the form 



E E 



1-2/ (vpSk^) 



kiieD2k±eDi ' 



less singular 



(38) 



(39) 



(40) 



where Di is a domain such that vpSk^ < T while D2 is the domain such that K:5k'^^ < T. The integration over D2 



thus gives a factor proportional to vT. The integral over Sk± can be transformed into an integral over energy in the 
same way as the perfect nesting case, with a constant density of states determined by the Fermi velocity. The domain 
delimited by Di and D2 is depicted in Fig. 3. 
Overall then, the final result will be that 



X 



(^)^^(Q,,c.)^TV2tanh(^) 



where the T-*-/^ prefactor comes from the k.. integration. A similar reasoning leads to 



X 



(1) 



(Qd,0)«Ti/2ln(^^^ +A 



(41) 



(42) 



which means that the regular temperature independent term represented here by A dominates at low temperature. 

we find that 



Repeating the same analogous arguments for ^ ^^^2^'^^ , 



^0^ 



2^1/2 



1 



2^2 2^3/2 

which implies from the definition of Fq Eqn.(15) and Eqs. (41-42) that 

P ^2rT.i/2 1 



(43) 



(44) 



the same result for Fq, within logarithmic corrections, as if we had perfect nesting. In three dimensions, the correction 
compared to perfect nesting is determined by the area spanned by dk\\, proportional to dk'^^ ~ T, hence we would have 

Results of numerical calculations shown in Fig. 4 confirm the power law temperature dependencies found above. 
At the actual filling Uc where the Fermi surface is tangent to the antiferromagnetic zone boundary, the power laws 
extend to low temperature. 
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2. Higher order and scaling 

Given that and Fq are now temperature dependent, we should check whether the smah q and smah u expansion 
of the denominator that lead to Eqn.(ll) is still valid. Normally, the expansion is of the form 



1^1^ L 



X"p(q> = TT^ . , „2^2 , ^„4^2 — 727T^' (45) 



with a a constant. Since the function falls on a scale ~ C ^ the higher order term oq"^^^ ~ ^ ^ can be neglected. 
However, our case is different. The coefficients of the expansion in powers of q are singular at T = 0. For example 

behave ^^^^g^ 



q=0 



^ where the comes from counting the powers of T associated with derivatives in 



Eqs.(29) to (31) and the T^/^ from the restriction to the k\\ integral as usual. Knowing the scaling of we can 



rewrite 7^ 



dqf 



~ ^0 7^ SO that we are left with 

q=0 



where a' is a constant. The susceptibility will preserve a scaling form as a function of q/T and uj/T if the scaling 
exponent is z = 1. Indeed, in that case ^ 1/T and since Fq 1/T, ~ 1/T^/^, the susceptibility becomes 

X'4(q, u;) ~ ^ Im fcv , _ .d^. (47) 

with h'^c'^d' constants. Each higher power of has an additional power of 1/T^ coming from the additional 
derivatives of the non-interacting susceptibility and the scaling form is preserved to all orders. For frequency, there 
are also higher order terms, (u/T)^ etc. Hence we have the general scaling form 

X'4(q,^) = -^5(|,^). (48) 

We check in the following section that this is consistent with the TPSC self-consistency condition. 



C. Quantum Critical behavior 

1. Correlation length, spin susceptibility and NMR relaxation rate 

The quantum critical behavior has been thoroughly studied in Ref. for the case where the Ornstein-Zernicke form 
is valid. This analysis does not apply here because of the more general form of the spin susceptibility obtained in the 
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FIG. 5: Log- log plot for the temperature scaling of the correlation length. The scaling is estimated from the width measured 
along one of the reciprocal lattice wave vectors (tt — qh) at various fractions of the maximum height. On (a), for U = 5.56, 
the critical doping corresponds to ric = 1.2007, where the Fermi surface is tangent to the antiferromagnetic zone boundary. 
The temperature scale is too small to detect possible logarithmic corrections. The results are consistent with z = 1. On (b), 
deviations from 1/T occur at low temperature because, for the chosen value [7 = 6, the calculation is at the critical doping 
1.20096, slightly away from the point where the Fermi surface is tangent to the antiferromagnetic zone boundary. Also shown, 
Cas obtained from the Ornstein-Zernicke form Eqn.(12) using Usp from the self-consistency relation Eqn.(7) 



previous section. Nevertheless, it is interesting to note that with the Ornstein-Zernicke form, one obtains 

r^^^- (49) 
J- 

for both perfect and pseudo nesting. Hence, simply taking into account the temperature dependence of Fq, we recover 
z = 1 scaling, namely 

e - ^- (50) 

Note that with a temperature independent value for Fq we recover the more usual result^^'""^ z = 2. 

The physics of the result for the correlation length is however quite different from the calculation with the Ornstein- 
Zernicke form. Indeed, in the latter case, it is the self-consistency relation Eqn.(7) that determines the temperature 
dependence of the correlation length. In the present case, we found that temperature dependence in the previous 
section without invoking the self-consistency relation. We will show in Sec. (IV D) that indeed in our case, the self- 
consistency relation leads to irrelevant corrections to the temperature dependence of the correlation length. 

The scaling of the correlation length can be obtained from numerical calculations by plotting, for example, the 
inverse of the width of the real part of the spin susceptibility at zero frequency measured at various fractions of the 
maximum value. For U = 5.56, the critical doping corresponds to ric = 1.2007 where the Fermi surface is tangent 
to the antiferromagnetic zone boundary. For that case. Fig. 5(a) shows that whether we measure the width at half- 
maximum or at some other fraction of the maximum, that width scales essentially as 1/T, with small deviations 
probably coming from the fact that we have not reached the asymptotic regime. We also show on this figure the 
correlation length ^as obtained from the Ornstein-Zernicke form Eqn.(12) using Usp from the self-consistency relation 
Eqn.(7). Deviations from the 1/T power law occur if we measure the width of the spin susceptibility too far in the 
tails, i.e. for a small fraction of the maximum (not shown). As demonstrated in Fig. 5(b), deviations from 1/T 
also occur at low temperature for a value of U (= 6 in the present example) where the critical point does not occur 
precisely when the Fermi surface is tangent to the antiferromagnetic zone boundary. 

To conclude this section, we show that one can easily obtain the temperature scaling for two more quantities. First, 
from the general form for the susceptibility Eqn.(45) used above, the static susceptibility at (7r,7r) scales as 

which can be checked directly numerically, or more simply deduced from the temperature dependent results for 
Eqn.(49), and ^g, Eqn.(43). 
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FIG. 6: On (a), log-log plot of the temperature dependence of the imaginary part of the self-energy at various color-coded 
points on the Fermi surface. The color code is in the inset. On (b), the local exponent is given as a function of angle and 
temperature. The points near the hot spot at — 7r/4 behave as T^^^ over the accessible temperature range. Calculations 
are done with U — 5.56, t' — —0.175 and t" — 0.05 at the quantum critical filling Uc — 1.2007. The lower figures are the 
corresponding results for U = 6, t' = —0.175 and = 0.05 at the quantum critical point n = 1.20096. Since in that case the 
Fermi surface is not exactly tangent to the antiferromagnetic zone boundary, the T^^^ behavior near ^ = 7r/4 is recovered only 
if the temperature is high enough that details of the Fermi surface cannot be resolved. The black lines on (b) and (d) are the 
curves defined by Tonset = VFSk±{0)/2 where Sk±{0) is the component of ki? — (7r/2,7r/2) parallel to VF(7r/2, 7r/2) at a given 
angle 0. 



Finally, the nuclear magnetic resonance (NMR) relaxation rate can be obtained from the two-dimensional 
version of the Moriya formula 

Tr' = T lim / 1^ |2 ^"^(q''^) ^2^ (52) 

uj^O J (jj 

where |Aq| is proportional to the hyperfine matrix element. Taking this as a constant and using the general scaling 
form Eqn.(48), a simple change of integration variable shows that 

= T lim / ^^^M1^2 ^ ^3/2^ (53) 

However, the integral over momenta q also contains contributions far from the peak in the susceptibility. There the 
susceptibility is essentially temperature independent. There is thus a Korringa contribution T^^ ~ T that is dominant 
at low temperature. 
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2. Self- energy 

To find the scaling of the self-energy, we use the scaling form of the susceptibility Eqn.(48) to rewrite the self-energy 
Eqn.(18) in the form 

Specializing to two dimensions and remembering the scaling of the Bose and Fermi functions with frequency and 
temperature, we change integration variables to x = ^ and y = ^ and we are left with 

S"^(k^,c.) = T3/25(^) (55) 

where S{^) is a scaling function. 




CO/t e/TT 

FIG. 7: On (a) and (c), log-log plot for the frequency dependence of —T^"^ (ki?,C(;) at various angles along the Fermi surface. 
The result of a power law fit is shown on (b) and (d). The dashed lines on (a) and (c) correspond to the fitted power laws. At 
the hot spot located at ^ = 7r/4, —T^"^ scales as cj^^^. The frequency range is small because of the low temperature saturation 
shown on the next figure. We have verified that the crossover from uj'^^'^ to the Fermi liquid regime uj'^ occurs on a broader 
angular scale when the temperature is higher, as expected from the results of Fig. 6. Calculations are done at T = 0.002, 
t' = -0.175 and t" = 0.05 at [/ = 5.56 and n = ric = 1.2007 for (a) and (b) and [/ = 6, n = 1.20096 for (c) and (d). 

The latter result can be checked numerically at = where we expect T^"^ (ki?,0) ~ 7^3/2 pjgg gj^g^^ Q^^id 6(c) 
displays the temperature dependent scattering rate for various angles 9 along the Fermi surface ioi U = 5.56 and 
n = Tie = 1.2007 and U = 6 and n = 1.20096, respectively. The line ^ = is horizontal in the Brillouin zone appearing 
in the inset. In Fig. 6(a), at the hot spot located at ^ = 7r/4, we recover the predicted result, T^/^. This is best 
illustrated in Fig. 6(b) by a plot of the local slope of the preceding log-log plot. As we move away from the hot spot, 
Fermi liquid behavior appears to be recovered. There are well known logarithmic corrections in two dimensions that 
may explain why we seem to deviate from exactly T^. One also notices that the T^/^ behavior occurs over a wider 
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FIG. 8: Log-log plot for the frequency dependence of —T^"^ at the hot spot for two different temperatures. The saturation at 
low frequency occurs at higher frequency when the temperature is higher. Calculations are done with [/ = 6, = —0.175 and 
t" — 0.05 at the quantum critical filling ric — 1.201 appropriate for electron-doped cuprates. 

range of angles when the temperature is high. This is easily understood from Figs. 6(b) and 6(d) that illustrates how 
temperature affects the domain where the pseudo-nesting occurs in the spin susceptibility. The solid black line in 
Figs. 6(b) and 6(d) is defined by Tonset = VFSk±{6)/2 where Sk±{6) is the component of k^? — (7r/2,7r/2) parallel to 
Vi?(7r/2, 7r/2) at a given angle 0. Figs. 6(c) and Figs. 6(d) are for U = 6,n = 1.20096. Since at this quantum critical 
point the Fermi surface is not tangent to the antiferromagnetic zone boundary, the T^/^ behavior occurs near = tt/A 
only at high enough temperature. At low temperature, deviations become apparent. 

When uo ^ the scaling form for the self-energy Eqn.(55) predicts S''^(kF,cj) uj'^' . However, at zero 
temperature, or when uj ^ the analytical approach taken above fails because the expansion of the spin susceptibility 
in uj/T and q/T is no longer justified and we cannot expect the latter result to be correct. Nevertheless, TPSC can be 
solved numerically. To set the stage for the next section where calculations are performed analytically directly at zero 
temperature, we show in Fig. 7 the result of the numerical calculations for a; T for two values of the interaction 
strength at a doping near Uc. On Figs. 7(a) and 7(b), U = 5.56 while U = 6 on Figs. 7(c) and 7(d). At the hot 
spot for U = 5.56, the scaling of the imaginary part of the self energy is very close to the expected result cj^/^. For 
[/ = 6, there is a larger discrepancy with the predicted scaling because at n = 1.20096 the Fermi surface does not 
touch (7r/2,7r/2) and thus the present theory does not apply anymore at low temperature. Away from the hot spot, 
Fermi liquid behavior is recovered. Again, logarithmic corrections are inaccessible from the numerical solution of the 
full TPSC equations because of the limited range of available frequencies: scaling is no-longer valid at frequency of 
the order of the Fermi energy, while at low frequency there is a saturation arising from the finite temperature. This 
saturation is illustrated in Fig. 8. We discuss analytically the T = regime in the following section where logarithmic 
corrections are found. The 1/a/T temperature dependence of the static (7r,7r) spin susceptibility obtained above will 
also be recovered. 



IV. FINITE FREQUENCY T = RESULTS FROM FIELD THEORY 
A. Lagrangian and scaling 

In this section we study the properties of fermionic excitations close to the hot spots within the field-theoretic 
framework of a spin-fermion model. This effective low-energy theory describes fermions with a parabolic dispersion 
(represented by fields, ip) interacting with SDW fluctuations (represented by a 0(3) vector field, 0). As shown in 
Fig. 1 there are four hot-spots on the Fermi surface which are connected by the SDW wave- vector Q = (tt, tt). Earlier 
studies^^'^^ of the spin-fermion model in the present context did not include the umklapp processes properly and we 
show in the following that a correct treatment of these terms modifies the results drastically. 

We start with the two patch (denoted by 5 = ±) model^^''^^ in the rotated {kx^ ky) coordinate system. The umklapp 
contributions will be discussed later. In order to simplify the notation, we have rescaled our coordinates to get rid of 
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FIG. 9: The polarization bubble for the two-patch theory in Eqn.56. The internal solid lines in the loop correspond to the free 
fermion propagators (different patches denoted by s = i) and the external wavy lines correspond to the boson 0". 

the Fermi velocity and curvature of the Fermi- surface. The corresponding Lagrangian takes the form 



f(V^T + f^- + fi^f. (56) 

Here we have promoted each fermion field to have N flavors (the flavor index is suppressed). The Yukawa- coupling, 
A, is chosen to be (9(1). As a result of this, the couplings of all the bosonic terms in the last line above are scaled 
by a factor of N as they will appear naturally upon integrating out the fermion fields. We don't include the kinetic 
energy of the boson, (9^-^)^, as this is an irrelevant term . 
The bare fermion propagator is given by, 

g°(fc)= ,fc +lk +fc2 ' ^ = (^-k). (57) 

The Fermi surfaces are located at = ky and kx = —ky for patch — and + respectively. From Fig. 1 we immediately 
observe that Q — (0, 27r) and also Q — (27r, 0) connect two more points in the BZ. However, it is convenient to fold 
back the points within the BZ, which effectively gives rise to two more patches. These can be described by rotating 
the original patches by 7r/2. Let us denote {ky^kx) by k. Then, the equations of these two additional fermi surfaces 
are given hy ky = k^.ky = —k^. 

Physically, these two scattering processes are very different since in the former case, the ^—fluctuation scatters 
fermions that disperse strongly in the direction transverse to the Fermi surface while in the latter case, they disperse 
strongly in the tangential direction. This will have interesting consequences in the behavior of the electron self energy 
as a function of the external frequency. 

The rest of this section is organized as follows: In section IV B, we compute the RPA contributions to the SDW 
propagator including both direct as well as umklapp processes. We then use the dressed bosonic propagator to evaluate 
the fermion self-energy in section IV C at leading order in 1/A^. 



B. RPA polarization 

At T = 0, the one loop polarization bubble (Fig. 9) for the two-patch theory is given by, 

U'^\q) = 2NS-»X^I ^[Gl{l)G^_{l + q) + Gl{l + q)G^_{l)], 
U{q) = iV[no(9r,q)+no(9r,-q)], 



(58) 
(59) 



where we are working with imaginary frequencies q^- and a, b denote the three SDW-polarizations. After performing 
the integrals, we obtain^^''^^ 



Ho(^r,q)-Ho(0,0) 



-Re 



:7r 



■i\qr 



El 



(60) 
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where, 



(61) 



In the RPA propagator obtained after bubble summation, the Ho (0,0) contribution determines the location of the 
quantum critical point. It is thus convenient to add and subtract this component to make the integrals convergent. 
From now on, we include no(0,0) in the definition of the bubble. 

To make connection with results of the previous sections, we also quote the results for the bubble at finite temper- 
ature. In this case. 



no(^) 



dz- 



1 



1 



where the contour C has to be chosen appropriately. Therefore, the above integral simplifies to. 



A2 



-Re 



TT 



1 



2x 



On integrating the above equation by parts, we obtain, 

A2 



where, 



no (9) = 
no (9) = 



J \ rj-1 5 rj-i 



167r2 



du- 
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cosh^ [{uj - Eq)/4T] 



167r2 



dy 



cosh'[{y-EjT)/4] 



2^2 



(62) 

(63) 

(64) 
(65) 

(66) 



as found in Ref. 

The RPA propagator for the SDW fluctuations is then given by. 



[no(<?r, q) + no((7r, -q) + ^oiQr, q) + ^oiQr, -q)] ' 



(67) 



where we have included the RPA contribution arising from all four hot spots on the Fermi surface and q = {qy^Qx)- 
The terms (no(^r5q) + ^o{Qt^ —^)) are not equal to {Ilo{qr^q) + Ilo{qr^ —q)) as was incorrectly assumed by the 
authors of Ref. ^. At the quantum critical filling r = 0, zero Matsubara frequency qr = but finite temperature, the 
term is negligible compared to the contribution of the bubbles IIq. Using the scaling form Eqn.(65) and keeping 
only terms linear in q in E^^ (Eqn. 61), one recovers the zero frequency limit of the previous result Eqn.(48) for the 
spin susceptibility. Naively doing the analytical continuation in frequency, the full scaling form would also follow. For 
the rest of the computations, we consider T = and carefully take into account logarithmic corrections that were 
beyond the reach of the previous calculation. 

Before we proceed to evaluate the fermionic self-energy, let us compute the forms of the real and imaginary parts 
of the retarded polarization bubble at T = 0, ni^(q, 1]), where i\qr\ ^ ^ -\- iO+. For the imaginary part, we obtain 
from Eqn. (60), 



Imn^(q,^]) 
Imn^(q,l]) 



A2 



2V^7r 
2V^7r 



Im 



^/E^ - - io+ + ^/E^ ^n^m 



0{n - E^)y/n-E^- o{-n - e^)^-q.-e^ 



(68) 
(69) 



which is chosen in a way such that f2Imnfi(q, f2) < 0. The real part can also be obtained from Eqn. (60) or from the 
Kramers-Kronig relation 



Renfi(q,i^) 



A2 



e{Eq - n) ^jE^-n + + ^q) ^j^ + e^ 



(70) 



These results agree with those of Refs. 22,37. The Vt <C limit calculated in Ref. also agrees with the above. 
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p-q 




FIG. 10: The electron self energy E+(p). The 0— propagator includes the one loop bubbles evaluated earlier. 

C. Electron self-energy 

We are interested in evaluating the electron self energy (Fig. 10) at T = 0, which at leading order in 1/A^ is 

given by, 



3A^ f (Pq 



(71) 



After analytic continuation iun uj-\- iO^ we obtain the following expression for the imaginary part of the retarded 
self-energy 

/72 
+ nsico - C^_q)]lmDfi(q, w - ^^_^), (72) 



where ^+ = + and = -kx + ky. 
At T = and k = we are left with 



ImE^ (0, Lo) = 3A2 I ^ [e(Cq) - e(Cq - w)] Imi?fl(q, to - C 



(Pq 
3A2 f , 



47^2/ "^^J/ / 2 '^9a;ImDi{(q,W-^_q) 



y d(7^ImDi^(q,cj-^_q) ^ (73) 

Note that = g^^^ + = for g^^^ — ^ — g^^. In terms of 11^^, this can be rewritten as, 

yRfn N 3-^^^ -Imn^^^(q,a;) 
lm^+(0,a;) - J dqy ^ Ken^^-^(q, c^))^ + (Imn^/(q, c^))^ ' ^^^^ 

where n^^(q, cj) is the total retarded RPA bubble including direct and umklapp terms and we have ignored the 
Qx Qy term in the denominator. 

For = -q^, we get = 3(7^/2, £;_q = E^^ = q^/2 - qy and £;_q = q^/2 + qy. Therefore, Re(Im)n^^ are 
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given by, 



Ren*?'(q,w) 



An 



^-2LO-ql0{-2uj-ql) + ^2lu - 6(210 - q^) 



+ ^q^-2qy-2uje{q^ - 2qy - 2oj) 



+ Jql + 2qy-2uj0{ql + 2qy-2Lo) 



Inin5?*(q,c^) 



47r 



^2uj-Zqie{2uj - 3ql) 



^q^-2qy + 2oj e{q*y - 2qy + 2oj) 
^J(lv + + 2u:e{ql + 2qy + 2w) 
^2u: + qie{2uj + ql) 



-3ql-2oje{-3ql-2u:) - J q^ - 2co Oiq^ - 2co) 



+ ^2uj + 2qy - q4 0(2uj + 2qy - q^) + ^2uj - 2qy - q^0{2uj - 2qy - q^) 
-^2qy -q^-2oj 0{2qy - q^ - 2cu) - ^ -2qy - q^ - 2cu e{-2qy - q^ - 2oj) 



(75) 



As a starting point, we can drop the Qy terms in the hmit of smah cj and retain only the qy terms in the q contributions 
that we take into account. This is a consistent way of handhng these terms, since if typical Qy ^ cu^ then Qy is smaller 
than Qy, so that it is justified to drop these terms. Since the integrand is an even function of g^, we integrate only 
over Qy > 0. Then, the expression for self energy reduces to. 



A 



B 



J — ( 



-Imn5?*(q,cj) 



^ (g2 + Ren5?*(q, uj)Y + (Imn5?*(q, i^)Y 



i dqyA{qy,uj) -\- / dqyB{qy,uj) -\- dqyC{qy,uj) -\- / dqyD{qy,uj) 

Jo Juj J ^2uj/2> JV2uj 



^2uj + 2qy + ^2uj - 2qy + J2uj - Sq^ + J2uj + ql 



A2 % + 47r 



1 2 



[y/2u;+2g^ + y/2u;-2g^ + y/2u;-3g2 + ^2u;+g2]: 
167r2 



2 I \/2g^+2u;+^2g^-2a;+^3gg+2u;+^2u;-gg 
A2 % + 47r 



[^2u;+2g^-y2g^-2a;+^2u;-3g2 + ^2u;+gg]- 
167r2 



c = 



^2a; + 2qy - ^J2qy - 2uj + ^2uj ^ q, 



j_ 2 , V2g^+2u;+y2g^-2u;+y/3gg+2a;+y3g2-2u;+y2u;-gg 
A2 % + 47r 



[y2u;+2g^ -y/2g^ -2u:+^2uj+ql 



167r2 



^2g^ + 2a; - ^2g^ -2uj- ^ql - 2a; + y^g^ _^ 2uj 



j_ 2 , v/2gy +2a;+y^2q^ -2a;+y3gg +2a;+y^3qg -2a; 
A2 % + 47r 



n 2 



[ v/2g^ +2a; - y^2g^ - 2u -^/^^+^ql +2a;] 2 
167r2 



(76) 



At small frequencies the dominant contribution to the imaginary part of the self-energy comes from the second integral 
between uj and ^2cj/3, which scales as ^ — a;^/^loga;. The contributions from the other regions scale as ~ o;^/^ and 
thus are negligible at low frequencies. The correct prefactor can be obtained by expanding the numerator and the 
denominator of integrand B for small frequencies uj and retaining only the largest terms, which gives 



B{qy,uj) 



47r'^V2^ 



(77) 
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FIG. 11: Left: ImSH(0, cj), calculated numerically from the second line of Eqn. (73), shown as a red solid line. The black 
dashed line is the asymptotic result Eqn. (78). For comparsion, the blue dash-dotted line indicates cj^^^. Right: different 
contributions to Im5]i?(0, cj). The red line shows the dominant contribution, denoted by B in the main text. Black dashed line: 
asymptotic result Eqn. (78). The other three curves show the sub- leading contributions A, C and D, which scale as ^ u^^^ . 



One then integrates over qy to obtain, 

We note here that the self-energy is less singular compared to earlier works where umklapp scattering was not taken 
into account, in which case the self-energy scales as — cjlogo;. ' . 

So far we haven't addressed an important issue: what happens if we include the renormalization of the boson- 
fermion vertex? In the two-patch theory originally considered by Altschuler et al. , the one-loop correction to this 
vertex was found to be logarithmically singular. However in Appendix A, we show that the full four-patch theory 
does not have this singularity. 

Based on our discussion so far, we see that the additional umklapp terms considered in our work play a very crucial 
role at the critical point. In the absence of these contributions, the self-energy was more singular than what we have 
found here. Moreover, the vertex correction was also found to be singular. However, here we have shown that the 
singular behavior is washed out when we include the additional scattering contributions. 



D. Effect of the self-consistency 

In the TPSC approach we imposed two-particle self-consistency in the form of a sum-rule that is similar to the 
spherical model. In the present field-theory approach, this amounts to imposing < 0^ >= 1 where the expectation 
value is taken with respect to the fermions and bosons. We have argued that the z — \ scaling does not come from 
the self-consistency condition. To confirm this result, in this subsection we obtain the scaling of the quartic term in 
the boson Lagrangian. 

Integrating out the fermions, the polarization operator Ho(g'r,q) appears in the quadratic term of the boson La- 
grangian. Emphasizing the scaling only, this term is symbolically written as 

.A 

y {d\duj)$^y^{uj,q). (79) 

This is the most relevant quadratic term. Integrating out the large wave number modes for q > A/s and rescaling q 
and uj such that q' = qs^ uo' = ujs with 5 > 1, returns the new cutoff A/s to its original value A. Invariance of the 
quadratic term written in terms of the prime variables then imposes that (j) = (p's"^ = (\)' 1^ . The effect of this tree 
level scaling on the quartic term is that 

u J {d\dLjf$^ ^ us-^s^"^ J {d^q'duj'f$'^. (80) 

This in turn means that u' = us~^s'^'^ = us~'^ scales to zero and is thus irrelevant. 
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V. SUMMARY 



We have argued that for bare interaction strengths U in the intermediate couphng range, commensurate SDW 
fluctuations at (tt, tt) and band parameters similar to those of electron-doped cuprates, the antiferromagnetic quantum 
critical point naturally occurs close to the fllling where the Fermi surface points joined by (tt, tt) are nearly tangent 
to each other. As long as the temperature or frequency are not too low, the limiting case of tangent Fermi surfaces 
describes the physics. In this pseudo-nesting situation, the Fermi liquid behavior breaks down. Quasiparticles still 
exist but the self-energy and spin susceptibility, for example, are different from those predicted by Fermi liquid theory. 

We considered this problem at zero temperature, or for frequencies larger than temperature, using a fleld-theoretical 
model of gapless collective bosonic modes (SDW fluctuations) interacting with fermions. The imaginary part of the 
retarded fermionic self-energy close to the hot spots scales as — o;^/^ logo;. This is less singular than earlier predictions 
of the form — uj log uj. The difference arises from the effects of umklapp terms that were not included in previous 
studies. 

At flnite temperature, we have used TPSC to study this problem and have obtained numerical results for the 
one-band Hubbard model with band parameters and interaction strength appropriate for electron-doped cuprates. 
Neglecting logarithmic corrections, we found analytically and numerically that the correlation length ^ scales like 
1/T, namely z = 1 instead of the naive z = 2. The static spin susceptibility x scales like 1/\/T, and the correction 
~ 2^3/2 Korringa NMR relaxation rate is subdominant. NMR experiments are diflicult in electron- 

doped cuprates. We also found that the imaginary part of the self-energy at the hot spot scales like T^/^. The 
latter result and the — cj^/^logo; frequency dependence of the self-energy should be experimentally veriflable with 
angle-resolved photoemission spectroscopy (ARPES) in electron-doped cuprates. Recent transport measurements in 
these compounds have found a T^/^ behavior of the resistivity above the quantum critical point at the end of the 
overdoped side of the superconducting dome. While there may be a relation with the above result if antiferromagnetic 
fluctuations disappear at the same time, one must also be careful not to equate scattering rate with resistivity because 
in a simple picture it is the inverse of the scattering rates that are averaged over the Fermi surface. This suggests 
that in the resistivity, Fermi liquid behavior of the cold spots should short-circuit the non-Fermi liquid behavior of 
the hot spots. 
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Appendix A: Vertex correction 

In this appendix, we compute the 1-loop correction to the Boson-Fermion vertex, which is deflned as, 

- {^P-AQ)^lAPmQ)) = r:,''i^^^_^liP,Q)i2^fSHQ-P-Q), (Al) 
where we are working again with imaginary frequencies. The expression for the diagram in Fig. 12 can be written as, 

/r]J 
-I—Gl{l+p + q)GUl+p)m- (A2) 

We now use the identity r^r^ = 5^'^ -{- ie^^^r^ twice to simplify the above expression. Then on deflning (^F^^, (p, q) = 
T^^;ST{p,q), we have. 



ST{p,q) 



N J (27r)3 
1 



-i{lr + Pr + Qt) + {Ix + + Qx) + {ly + P^y + Qy)^ 

1 



-i{lr + Pr) - (Ix + Px) + {ly + PyY 



p + r + [HoO) + no(-i) + no(i) + no(-i)] 



(A3) 
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FIG. 12: The 1-loop contribution to the Boson- Fermion vertex. 



Let us now evaluate this for zero external momenta q = p = at the critical point and check for singularities. The 
expression reduces to, 



. _A3 f dlrdH 

N J (27r)3 



1 



ilj- ~\~ Ix ^'y 



1 



irj. -\- ly 



P + Uo{l, Ir) + no(-l, Ir) + no(l, Ir) + Uo{-l, Ir) i 



(A4) 



The above integrand has a very complex structure. Let us therefore analyze the (non-) singular nature of this diagram 
by power counting. One needs to be careful as the ^—propagator has many combinations of powers of the internal 
momenta. We begin by rescaling the variables as, 



/ — / / — /^/ / 



2,1 



(A5) 



Eqn. A4 then takes the form, 

5V 



N J (27r) 



dl 



-il'^ + + l 



1 



-ilL 



I'+l 



^ {i%y2 ± 1) + Vii%y2 ± 1)2 +7^ 



(A6) 

(A7) 
(A8) 



Ah we need to do now is to check whether this expression (which is so far exact) is singular in the IR and UV. In the 
IR, we can ignore the momentum dependence of the fermionic Green's functions compared to 1 in the denominator. 
Moreover, f{l'^,l'^) « 1 and g+{l,l'^,l'^) « l,g-{l,l'^,l'^) « in this small momentum limit. Therefore, the above 
expression reduces to. 



,5r 



A3 

' N 



r dl'^dl'^ dl 

Jo 



P + {2l + Vl)/2Tr' 



(A9) 



where e is a small cutoff. But the above expression is convergent, so that there are no IR singularities. 

Let us now check for UV singularities. We proceed by introducing a characteristic lower cutoff A which is large, 
but finite, such that the integration runs from A to co. Then, in the limit of these large momenta, we have. 



5T 



A3 r°° dl'^dl'^dl 



1 



-il' + 1' 
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-ill - I' 
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p + iij^w+w+i'.+'^mi'^T^ 



(AlO) 
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where we have ignored the contribution from f{—l'^^l'^) compared to f{l'^^l'^) and I'^l'^ compared to l^l'^- We first 
want to do the integral over l'^ and l'^. It is more convenient to change to polar coordinates, (/^, l'^) (r, 6). However, 
we estimate 1'^. ^ r and eliminate the 9 dependence, which does not give rise to any singularities. Then, 



Both the r and / integrals can be performed easily to verify that 5T is convergent in the UV. Therefore, the 1-loop 
vertex correction is non-singular both in the UV and in the IR. 
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